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Abstract 

We have recently solved the inverse scattering problem for one pa- 
rameter families of vector fields, and used this result to construct the 
formal solution of the Cauchy problem for a class of integrable nonlin- 
ear partial differential equations connected with the commutation of 
multidimensional vector fields, like the heavenly equation of Pleban- 
ski, the dispersionless Kadomtsev - Petviashvili (dKP) equation and 
the two-dimensional dispersionless Toda (2ddT) equation, as well as 
with the commutation of one dimensional vector fields, like the Pavlov 
equation. We also showed that the associated Riemann-Hilbert inverse 
problems are powerfull tools to establish if the solutions of the Cauchy 
problem break at finite time, to construct their longtime behaviour 
and characterize classes of implicit solutions. In this paper, using the 
above theory, we concentrate on the heavenly and Pavlov equations, 
i) establishing that their localized solutions evolve without breaking, 
unlike the cases of dKP and 2ddT; ii) constructing the longtime be- 
haviour of the solutions of their Cauchy problems; iii) characterizing 
a distinguished class of implicit solutions of the heavenly equation. 
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1 Introduction 



It was observed long ago pQ that the commutation of multidimensional vector 
fields can generate integrable nonlinear partial differential equations (PDEs) 
in arbitrary dimensions. Some of these equations are dispersionless (or quasi- 
classical) limits of integrable PDEs, having the dispersionless Kadomtsev - 
Petviashvili (dKP) equation [2], [3] as universal prototype example; they arise 
in various problems of Mathematical Physics and are intensively studied in 
the recent literature (see, f.i., [I] - [2 1 J ) • In particular, an elegant integration 
scheme applicable, in general, to nonlinear PDEs associated with Hamil- 
tonian vector fields, was presented in [B] and a nonlinear d - dressing was 
developed in [12]. Special classes of nontrivial solutions were also derived 
(see, f.i., HI], M)- 

Distinguished examples of PDEs arising as the commutation conditions 
[Li(A), L 2 (X)} = of pairs of one parameter families of vector fields, being 
A G C the spectral parameter, are the following. 
1. The vector nonlinear PDE in N + 4 dimensions 



u tlZ2 - u t2Z1 + (u zi ■ v*) u Z2 - (u Z2 ■ V*) U Z1 = 0, (1) 

where U(ti,t 2 , z 1 ,z 2 , x) E R N , x = (x 1 ,..-,^) G R N and = (d Xl , .., d XN ), 
associated with the following pair of (N + 1) dimensional vector fields 

U = d u + \d Zi + U Zi ■ Var, i = 1,2. (2) 
2. Its dimensional reduction, for N = 2 1221: 



U tx - U zy + (U v ■ V^j U x - (U x ■ Vaf) U y = 0, 
UeR 2 , x= (x,y), Vg=(d x ,d y 



(3) 



obtained renaming the independent variables as follows: t\ = z, t 2 = t, x\ 
x, x 2 = y, associated with the two-dimensional vector fields: 

Li = d z + \d x + U x ■ Vj?, 
L 2 = d t + \d y + Uy ■ V s . 



(4) 



3. The divergenceless reduction ■ U = of (J3J), the celebrated second 
heavenly equation of Plebanski [23] : 

Oir - ■ JK m - 0% 0. 6 = 6(x,y,z,t) eR, x,y,z,teR, (5) 
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describing self-dual vacuum solutions of the Einstein equations, associated 
with the following pair of Hamiltonian two-dimensional vector fields 



(6) 



(7) 



(8) 



L\ = d z + Xd x + 6 xy d x — 6 xx d y , 

L2 = 9 t + Xdy + Oyyd X ~ Q X y3 y . 

4. The following system of two nonlinear PDEs in 2 + 1 dimensions \2 L 

1l x t ~\~ Uyy ~t~ (tlU x ^ x ~\~ V x U x y VyU xx 0, 
V X t Vyy ~\~ UV XX ~\- V x V x y VyV xx 0, 

arising from the commutation of the two-dimensional vector fields 

L x = d y + (A 4- v x )d x - u x d x , 
L 2 = d t + (A 2 + Xv x + u- v y )d x + {-\u x + u y )d\, 

and describing a general integrable Einstein- Weyl metric |25j. 

5. The v = reduction of ([7]), the celebrated dKP equation 

(u t + uu x ) x + u yy = 0, u — u(x, y, t) G R, x,y,teR, (9) 

(the x-dispersionless limit of the celebrated Kadomtsev-Petviashvili equa- 
tion [26J), associated with the following pair of Hamiltonian two-dimensional 
vector fields [51 E] : 

L 1 = d y + Xd x - u x d x , / 1Q x 
L 2 = d t + (A 2 + u)d x + ( Xu x + Uy)d\, 

describing the evolution of small amplitude, nearly one-dimensional waves in 
shallow water [27] near the shore (when the x-dispersion can be neglected), 
as well as unsteady motion in transonic flow [2] and nonlinear acoustics of 
confined beams [3]- 

6. The u = reduction of (J7J), the Pavlov equation [T5] 

Vxt + Vyy = VyV xx ~ V x V x y, V = V (x , ?/ , t ) G R, X, |/, t G R, (ll) 

associated with the non-Hamiltonian one- dimensional vector fields 

Li = d y + (A + v x )d x , , . 

L 2 = d t + (X 2 + Xv x -v y )d x . 1 j 



7. The two-dimensional dispersionless Toda (2ddT) equation [2E1 [2H] : 

fe. = = 0(Ci,C 2 ,t) (13) 

(or (^Cif2 = (, eV )tti V 9 = 0t)> associated with the pair of Hamiltonian vector 
fields 0: 

^ = % + Ae^ft + {-\{e% + %±) A«9 A , 

L 2 = «9 ?2 + A-^ft + (A- 1 ^), - %i) A«9 A , (M) 

describing integrable heavens [301 EI] an d Einstein - Weyl geometries |32j, 
[53] , [Mj; whose string equations solutions [S] are relevant in the ideal Hele- 
Shaw problem [35]-[39|. 

The Inverse Spectral Transform (1ST) for 1-parameter families of multi- 
dimensional vector fields, developed in [22] (see also [10]), has allowed one to 
construct the formal solution of the Cauchy problem for the nonlinear PDEs 
(|3D and (J5D in |22J, for equations © and © in [21], for equation ffTT]) in [4T] 
and for the wave form (e^')t = 4> xx + 4> yy of equation ( fl3l) in [42J. This 1ST, 
introducing interesting novelties with respect to the classical 1ST for soliton 
equations [SI [27], turns out to be, together with its associated nonlinear Rie- 
mann - Hilbert (RH) dressing, an efficient tool to study several properties of 
the solution space of the PDE under consideration: i) the characterization of 
a distinguished class of spectral data for which the associated nonlinear RH 
problem is linearized and solved, corresponding to a class of implicit solutions 
of the PDE (for the dKP and 2ddT equations respectively in [45J and in [42J , 
and for the Dunajski generalization [16] of the heavenly equation in |44J); ii) 
the construction of the longtime behaviour of the solutions of the Cauchy 
problem (for the dKP and 2ddT equations respectively in [15] and [12]); iii) 
the possibility to establish whether or not the lack of dispersive terms in the 
nonlinear PDE causes the breaking of localized initial profiles (for the dKP 
and 2ddT equations respectively in [15] and in [12]) and, if yes, to investigate 
in a surprisingly explicit way the analytic aspects of such a wave breaking, 
as it was done for the dKP equation in [15] . 

In this paper we first develop, in §2, the RH dressing schemes for equa- 
tions (jSJ) and ffTTj) . Then we use these two schemes 

i) to show that, in contrast to the cases of the dKP and 2ddT equations, 
localized solutions of the PDEs (jSJ) and (fTTj) do not break (in §2); 

ii) to construct the longtime behaviour of the solutions of their Cauchy prob- 
lems (in §3); 
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iii) to characterize a distinguished class of spectral data for which the non- 
linear RH problem of equation <^ is linearized and solved, corresponding to 
a distinguished class of implicit solutions of (jSJ) (in §4). 



2 Riemann - Hilbert dressing 

In this section we present the RH dressing schemes for equations and ffTTj) . 
and show that localized solutions of these two models do not break. Both 
schemes can be extracted from the inverse spectral transforms presented in 
[22] and in [41] . In solving the Cauchy problems for equations (J5j) and ffTTj) . 
the RH data of these dressing schemes are connected to the initial conditions 
via the direct spectral transforms presented in [22] and in [41] . 

1. RH dressing for (EI). Consider the vector nonlinear RH problem on the 
real line: 

7r+(A) = tT(A) +R(tT(\),\), AgR, (15) 

where 7? + (A), 7?~(A) G C 2 are 2-dimensional vector functions analytic in the 
upper and lower halves of the complex A plane, normalized as follows 

7 f ± (A) = v(\;x,y,z,t) + 0(A~ 1 ), |A| » 1, 

x - \z \ (16) 



v(\;x,y,z,t) -j y _ Xt 

and the spectral data R((, A) = (i?i(Ci, C2? A), i?2(Ci> C2, A)) G C 2 , defined for 
C G C 2 , A G K, satisfy the following properties: 



TZ(JZ((, A), A) = C, VC G C 2 , reality constraint, 
{7?-i, 7^2},? = 1, heavenly constraint, 

where TZ((, A) := C + -R(C> A) and {•, is the usual Poisson bracket with 

respect to the variables ( = ((1,(2)- Then, assuming uniqueness of the 
solution of such a RH problem and of its linearized version, it follows that 
are solutions of the linear problems Li^if = 0, where L 12 are defined in 
(MI), and 

_^ J =F(x,y,z,t) GK 2 (18) 
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is solution of the heavenly equation (151) , where 

/d A 
— R\iri(\; x, y, z, t), tt^(A; x, y, z, t),Xj. (19) 

K 

As a consequence of equations Li^^ = 0, it follows that, for |A| >> 1: 

nf = x - \z - 9 y X~ l + 9 t \- 2 + 0(A- 3 ), 

7r± = y - At + ^A" 1 - # 2 A~ 2 + 0(A- 3 ). 1 ' 

Proof. We apply the operators in (jlj) to the RH problem ffl5l) . where 

(»:«))• < 21 > 

obtaining the vector equation Ljn + = J(ir~, X)Lj-ix~, j = 1,2, where J is 
the Jacobian matrix of the transformation f[T51) : Jfc/(C, A) = <97£fc(C, X)/dQ. 
Since, due to (I2~T|) . — > as A — > oo, by uniqueness we infer that ^ 
are eigenfunctions of the vector fields Lf LjTT^ = 0. Consequently, U is 
solution of equation ([3]). In addition, if the reality constraint (fTTK ) is im- 
posed, then, by uniqueness, it follows, from (TIB]) , that 7r+ = 7T~, implying 
the reality of U. At last, if the heavenly constraint (117b ) is imposed, then 
{nf,^} = {^1,^2}- Since {7rf,7r^} ^ 1 as |A| — > 00, the analyticity of 
the eingenf unctions implies that {n^ , n^} = {7rf , n 2 } = 1. Applying Li ]2 to 
these two equations, one infers that V# • ?7 = 0, implying the existence of a 
potential 9 such that U = (9 y , —9 X ), and implying also that the two vector 
fields L12 are Hamiltonian, with Hamiltonians (H 1 ,H 2 ) = Vg9, reducing to 
the vector fields (jHJ). Then the system ([3]) reduces to the heavenly equa- 
tion ([5]) and equation ( fl8|) follows from ( I2T1) and from the integral equations 
characterizing the RH problem. □ 

2. RH dressing for (II lft . Consider the scalar nonlinear RH problem on 
the real line: 

7r+(A) = 7r~(A) + J R(vr~(A),A), A e R, (22) 

where 7r + (A), 7t~(A) G C are scalar functions analytic in the upper and lower 
halves of the complex A plane, normalized as follows 

n ± (\;x,y,t) = u( y \;x,y,t)+0(\- 1 ), |A| » 1, _ 
z/(A; x, y, t) = — XH — Xy + x, 
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and the spectral datum R((, A) G C satisfies the following reality constraint: 



n(n(C,X),X) = C, VCeC, XeR, (24) 

where A) := C + -R(C>^)- Then, assuming uniqueness of the solution 
of such a RH problem and of its linearized version, it follows that 7r =1= are 
solutions of the linear problems L\$k = 0, where are defined in (fl~2l . 
and 

v — F (x,y,t) Gl (25) 
is solution of equation (ITT1) . where 

—R(7T-{X;x,y,t),Xy (26) 

R 

Proof. We apply the operators in l fi"2l to the RH problem ( l22l) . where 
u(x, y,t) = - lim A (vr ± + A 2 t + Xy - x) , (27) 

A— >oo 

obtaining I/j7r + = J(n~, X)Lj7r~, j = 1,2, where J((, A) = d7Z((,X)/d(. 
Since Livr^ 1 — >• as A — > oo, it follows that, by uniqueness, 7r =1= are eigenfunc- 
tions of L\. L\tx ± = 0. Consequently, we have the asymptotics 

7i ± = -X 2 t-Xy + x-l + ! + 0(A- 3 ), |A|»1. (2g) 

/3 x =Vy+ 

implying that also Z/27r ± — > as A — > oo. Therefore, again by uniqueness, 
^i^^ = 0, f satisfies equation ffTTl) and equation (1231) follows from (|2T|) and 
from the integral equations characterizing the RH problem. If, in addition, 
the reality constraint (|24j) is imposed, then, by uniqueness, 7r + = n~ and, 
consequently, t)6l □ 

We recall that, in the dressing schemes for the dKP and 2ddT equations, 
the spectral mechanism responsible for breaking is the presence, in the nor- 
malizations of the solutions of their RH problems, of the unknown fields, so 
that the inverse formulae define the solutions of the PDEs only implicitely 
[13,112]. Since the normalizations f|T6l) . fl23l) depend only on the space-time 
independent variables, such a breaking mechanism is absent here, and we 
expect that localized solutions of ([5]) and fTTT]) do not break during their 
evolution. 
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3 Longtime behaviour of the solutions 

The 1ST provides an effective tool to construct the longtime behaviour of 
the solutions of integrable PDEs. In this section we investigate the longtime 
behaviour of the solutions of the nonlinear RH problems (fl5l) . ([221 . and, 
consequently, of the solutions of the Cauchy problems for equations ([5]) and 

CCD- 

Asymptotics of equation §5$). We study the longtime t » 1 regime in 
the space regions 

x = x + fit, y = v 2 t, z = v 3 t, x, v i, t> 2 , v 3 = 0(1), t> 2 , t>3 7^ 0, t » 1. (29) 

The system of nonlinear integral equations characterizing the solutions of 
the RH problem (fi~5|) is conveniently written in the form: 

= h k A^feo)^ + (*i - A't*)t + 0i (A'), 
(*, 2 -A')t + 2 (A'),A'), j = 1,2, 

where 

0i(A) = vrr(A) - (a: - Xz), 2 (A) = vr^A) - (y - At). (31) 

The fast decay, in t, of fy, j = 1,2, due to the two different linear growths 
in t of the first two arguments of Ri, i? 2 , is partially contrasted if V\ = V2V3, 
i.e., on the saddle surface: 

x = x + y (vi = v 2 v 3 ). (32) 

Indeed, on such two-dimensional manifold, the first two arguments of R±, R2 
in (|30|) exhibit the same linear growth in t: 



<MA) = h k A^ffco)^ + HV2 - X')t + fc(A') 
(f 2 -A')t + 2 (A'),A'), j = 1,2. 



Since, in this case, the main contribution to the integral occurs when A' ~ v 2 , 
we make the change of variable A' = i> 2 + n' /t, obtaining: 

(*) = 2kJ ^A-(^-io) ^ (* - ^ + frfa + t)» , x 

k v (34) 

-^ + 02^2 + ^),^), .7 = 1,2. 
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If |A - v 2 \ >> t 1 , equation ([MD implies that (j)j(X) = O^T 1 ), j = 1,2: 
<M A ) ~ - art(A-(^+io))t I d ^' R J ( £ ~ + 0i (u 2 + £), 

; » v (35) 

V + 02^2 + ^),^), J = 1,2, 

while, for X — v 2 = fi/t, fi = 0(1), we have that 0j(A) = 0(1), j = 1, 2: 

fcte + f ) ~ h I w^m R 3 (s - «w' + Mv2 + £), 

k V (36) 

-^ + 02^2 + ^), V 2 ), J = 1,2. 



Therefore it is not possible to neglect, in the integral equations (Ell), 01,2 in 
the arguments of and these integral equations remain nonlinear also in 
the longtime regime. Summarizing, due to equations (1361) . (ITBl) and (fl9l) . the 
longtime behaviour of the solutions of the Cauchy problem for the heavenly 
equation on the saddle (132!) reads: 

( ~l V ) = i^oo (x - f , f , f ) + oir 1 ), x - f = 0(1), t » 1, (37) 



where 

F 0O (x,v 2 ,v 3 ) 



—±r J dfiR[x - v 3 /i + ai(ji; x,v 2 ,v 3 ), -/x + a 2 (/j; x,v 2 ,v 3 ),v 2 ) . 

TO V / 



(38) 



and <%,(/!; x, v 2 , v 3 ), j = 1, 2 are the solutions of the following nonlinear inte- 
gral equations 

R 7 

(39) 

Outside the saddle, the solution decays faster than 1/t. 

We remark that the longtime behaviour ( 1371) is formally the same as that 
of solutions of the linearized heavenly equation 9t x — zy = 0; the nonlinearity 
manifests only in the fact that is obtained, from the scattering data, 
through the solution of the nonlinear integral equations ([39]). These equations 
characterize the new RH problem: 

l + (/i) = A-(n) + R(A-(fj,),v 2 ), /iGK, (40) 



9 



where A ± (^l; x, v 2 , v%) G C 2 are 2-dimensional vector functions analytic in the 
upper and lower halves of the complex /x plane, normalized as follows: 

^)=( £ -^)+ogo, (4i) 

R(C, A) are the heavenly spectral data and a(/x; x, v 2 , v 3 ) = A~(fj,; x, v 2 , v 3 ) — 
(x - v 3 fi, -/x) T . 

The above analysis extends, with minor modifications, to equation ffTTj) . 
for which we give just the results. 

Asymptotics of equation ( 1111) . Let t » 1 and 

x — x + Vit, y = v 2 t, x,vi,v 2 = 0(1), v 2 ^ 0, t » 1. (42) 
On the parabola 

* = (^ = -f)' ( 43 ) 



the analogue of the heavenly saddle (1321) . the longtime behaviour of the 
solutions of equation (TTTT) is given by 



- = 7i F -{ x+V i^)+o(^), (44) 

where 

i^oo {x, «a) = 2^7 / - /i 2 + a(/x; x, v 2 ), -f ) , (45) 

-R(C, A) is the spectral datum of equation (TTTT) , and a(/x; x, u 2 ) is the unique 
solution of the nonlinear integral equation 

= 35? / ^fco) * (* + ^ + «(a0, f ) , (46) 

characterizing the solution of the following scalar nonlinear Riemann problem 
on the real axis: 

A+( f i)=A-(ji)+R(A-(ji),-f), /2 eR, 
A ± (fx;l;,f)=x-fx 2 + 0(fi- 1 ), |//| » 1, 1 ' 

where a(fi;x, y) = A~(/i) — (x — /U 2 ). Outside of the parabola (S3]), the 
solution decays faster than \j\ft. 
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4 Implicit solutions of the heavenly equation 

In this section we construct a class of explicit solutions of the vector nonlinear 
RH problem (fTB1) - (iT7j) and, correspondingly, a class of implicit solutions of 
the heavenly equation (jHJ), parametrized by an arbitrary real function of two 
variables. 

Suppose that the RH data in (ITS"]) are given by: 

RjiCu C2, A) = C^H^/CC^a) _ 1); J = 12) (4g) 

in terms of the new real spectral function /(£, A) depending on tt^" and 7r^~ 
only through their product. Then the RH problem (Tl5i) becomes 



7r+ = 7rre^K^^, A G R, 
and the following properties hold. 

i) The reality and the heavenly constraints (I17p are satisfied. 

ii) 71^712 = ^1^2 ■ Consequently, {^t^t) is just a polynomial in A: 

71^2 = ^1^2 = Z ^ ~ fa + U z )^ + xy — z9 x + t6 y = w(X), (50) 

and the vector nonlinear RH problem (|15p decouples into two linear scalar 
RH problems: 

7r + = 7 r 1 - e i/WA),A) j 



(51) 



iii) Since, from (|5"Tj) . 



3 3 

where 



7r +e«(-) i ^W = w - c «(-) i ^W > j = 1,2, (52) 



= ^ / - — ^— — /(w(A'),A'), (53) 
2m J R A — (A ± zO) 

also 

7r + e i(-y/+(A) ) j = i ; 2 are polynomials in A. We expand them in powers 
of A, for |A| >> 1, and introduce the notation 

< X n f >= -L / A"/ HA), X)dX, neN. (54) 
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From the positive power expansions, it follows that 

n + e -if+(x) = n - e -if-(\) =x _ z \_ z< f >; 



7T. 



+ e i/+(A) = n - e if~W =y-t\ + t< f >, 



(55) 



implying the following explicit solution of the RH problem fl5Tl) : 

^ = (x-Xz-z<f>)e^\ 

n± = (y-\t + t<f>)e-^ x \ [ ' 

From the coefficients of the A -1 term, which must be zero, one finally obtains 
the two conditions 

B y = x < f > -z (< Xf > +| < / > 2 ) , 

9 x = y<f>-t{<Xf>-\<f>*). v f > 

Since < f > and < Xf > depend, through w, on 9 X and 9 y (see fl50l) ). it 
follows that (157p is a system of algebraic equations defining implicitely the 
solution (9 x ,9 y ) of the heavenly equation. 

In addition, equating to zero the coefficients of the A -2 terms (using (1201)). 
one obtains: 

9 t = 9 y < f > -x (< Xf > +\ < f > 2 ) + z (< A 2 / > + 

< / >< Xf > +| < / > 3 ) , 

9 Z = -9 X < f > -y (< Xf > -\ < f > 2 ) + t (< A 2 / > - [ ™> 

< f >< Xf > +| < / > 3 ) , 

and manipulating equations fl58l) and fl57|) . one shows that the solutions of 
(j3J) generated by the spectral data (j48p satisfy the following linear PDE 

td t - z9 z + y9 y - x9 x = 0. (59) 



Substituting in (J5J) the expression of 9 t given in fl59|) . one obtains the lower 
dimensional nonlinear constraint 

(x9 x ) x - y9 xy + z9 xz - t9 yz + t(9 xx 9 yy - 9 2 xy ) = (60) 

satisfied by the solutions of (jSJ) generated by the spectral data (jUJ). 
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